Rotational dynamics of a superhelix towed in a Stokes fluid 
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Motivated by the intriguing motility of spirochetes (helically-shaped bacteria that screw through 
viscous fluids due to the action of internal periplasmic flagella), we examine the fundamental fluid 
dynamics of superhelices translating and rotating in a Stokes fluid. A superhelical structure may 
be thought of as a helix whose axial centerline is not straight, but also a helix. We examine 
the particular case where these two superimposed helices have different handedness, and employ 
a combination of experimental, analytic, and computational methods to determine the rotational 
velocity of superhelical bodies being towed through a very viscous fluid. We find that the direction 
and rate of the rotation of the body is a result of competition between the two superimposed 
helices; for small axial helix amplitude, the body dynamics is controlled by the short-pitched helix, 
while there is a cross-over at larger amplitude to control by the axial helix. We find far better, 
and excellent, agreement of our experimental results with numerical computations based upon the 
method of Regularized Stokeslets than upon the predictions of classical resistive force theory. 

PACS numbers: 



INTRODUCTION 

The study of swimming micro-organisms, including 
bacteria, has long been of scientific interest [j], 2, El- 
Bacteria swim by the action of rotating, helical flagella 
driven by reversible rotary motors embedded in the cell 
wall Q ■ Typically, these flagella visibly emanate from the 
cell body. The external flagella of rod-shaped bacteria, 
such as E. coli, form a coherent helical bundle when rotat- 
ing counter-clockwise, causing forward swimming. When 
these flagella rotate in the opposite direction, the flagel- 
lar bundle unravels, causing the cell to tumble. This run 
and tumble mechanism allows a bacterium to swim up a 
chemo-attractant gradient as it senses temporal changes 
in concentration [J, [Bj]. Many studies have focused on 
the fundamental fluid mechanics surrounding this loco- 
motion affected by a simple helical flagellum attached to 
and extruded from the cell body [H, Q . Recently, there 
have been additional studies that investigate the hydro- 
dynamics of flagellar bundling @, [1| . 

In contrast, swimming bacteria with more complicated 
body-flagella arrangements are less studied. Spirochetes 
are such a group of bacteria. They have a helically- 
shaped cell body [l(J, and although they also swim 
due to the action of rotating flagella, these do not visibly 
project outward from their cell body. Instead, the cell 
body is surrounded by an outer sheath, and it is within 
this periplasmic space that rotation of periplasmic flag- 
ella (PFs) occurs. These helical periplasmic flagella em- 
anate from each end of the cell body, but rather than 
extend outwards, they wrap back around the helical cell 
body. In the case of Leptospiracaeae, there are two PFs, 
one emerging from each end of the cell body, that do not 
overlap in the center of the cell. Rotation of each flag- 



ella is achieved by a rotary motor embedded in the cell 
body. The shapes of both ends of the helical cell body 
are then determined by the intrinsic helical structure of 
the periplasmic flagella, as well as their direction of ro- 
tation. During forward swimming, L. Mini exhibit an 
anterior region that is superhelical, due to this interplay 
of helical cell body and helical flagellum [Io| . In fact, the 
handedness of these two helical structures are opposite, 
with the flagellum (axial helix) exhibiting a much larger 
pitch than the cell body helix. 

The overall swimming dynamics of spirochetes involves 
non-steady coupling of the complex geometry of the cell 
body, the flexible outer sheath, and the counter-rotation 
of the cell body with the internal flagella. However, a 
natural question is how the effectiveness of spirochete lo- 
comotion depends upon the detailed superhelical geome- 
try of the anterior region of the bacterium. With this as 
motivation, we present here a careful study of the funda- 
mental fluid mechanics of superhelical bodies translating 
and rotating through highly viscous fluids. We extend 
the classical analytic and experimental results of Pur- 



eed [6] and the numerical results of Cortez et al. [11 1 
performed for regular helices. In addition, we offer co- 
ordinated laboratory and computational experiments as 
validation of the method of Regularized Stokeslets for 
zero Reynolds number flow coupled with an immersed, 
geometrically complex body. This method uses modified 
expressions for the Stokcslct in which the singularity has 
been mollified. The regularized expression is derived as 
the exact solution to the Stokes equations consistent with 
forces given by regularized delta functions. 

We focus on a typical body that is a short-pitched helix 
whose axis is itself shaped as a helix of larger pitch and 
opposite handedness. In the following sections, we de- 
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scribe the experimental set-up as well as the construction 
of these superhelical bodies. We experimentally measure 
the rotational velocities of the bodies as they are towed 
with a constant translational velocity through a very vis- 
cous fluid. Note that rotational and translational veloc- 
ities should be proportional, with the constant of pro- 
portionality (resistance coefficient) dependent upon the 
body geometry. The rotational velocities corresponding 
to translational velocities are also predicted analytically 
using resistive force theory, as well as using the method 
of Regularized Stokeslets [llj, |l2j. We find compatible 
behavior between experiments and the resistive force the- 
ory, but excellent quantitative agreement between exper- 
iments and the method of Regularized Stokeslets. 



SUPERHELIX CONSTRUCTION 

A superhelix is formed from a copper wire chosen to 
be sufficiently malleable to deform into a desired shape, 
but rigid enough not to deform as it moves through the 
viscous fluid. The superhelix is made in two steps (see 
Fig.[ljb)). First, a copper wire of diameter 0.55 mm is 
wound tightly in a clockwise direction up a rod of di- 
ameter 3.15 mm, forming a tight coil. After removing 
the coil from the rod, we stretch it out into a smaller 
radius, larger pitch helix, simply by pulling the ends of 
the coil away from each other. The axial helix is made 
in the same manner, but we use lead wire of a thicker 
diameter (3.15 mm) and a larger rod (4.7 mm diameter). 
The most important difference between the two helices 
is handedness; the axial helix is wound counter-clockwise 
up the rod, whereas the small helix is wound clockwise. 
Once the parameters of the small and axial helices are 
measured, the axial helix is threaded through the small 
helix, forming a superhelix, i.e. the small helix is placed 
back on a rod that has been distorted into a helical shape. 
The last step is to remove the axial helix. This is done 
by simply rotating the axial helix while keeping the su- 
perhelix fixed. 

The defining geometric parameters of the superhelix 
are the radius r and pitch p of the small helix, and the 
radius R and pitch P of the axial helix (see Fig.[T|). For 
R = or P = oo, the superhelix reduces to a regular 
helix. In our experiments, two different sets of small he- 
lices are used. The corresponding geometric parameters 
of these small helices are the pitch (5.58 ±0.25 mm; Set I 
and 5.04 ±0.36 mm; Set II) and radius (1.91 ±0.14 mm; 
Set I and 1.75 ± 0.21 mm; Set II). The small (less than 
12%) variations of pitch and radius are presumably due 
to mechanical relaxation of material when it is pulled off 
the axial helix. Seven different axial helices are prepared 
from the same initial coil (see Fig. [2|). 

We now construct a mathematical representation of 
the superhelix. The coordinates of an axial helix are X = 
(Rcos(Kz), Rsm(Kz), z), where K = The distance 
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FIG. 1: (a) Schematic of experimental setup. A motor pulls 
a rigid body through silicon oil, a highly viscous Stokes fluid 
(y = 10,000 cSt) (b) Procedure for making superhelix. 



measured along this helix is linearly proportional to the 
axial distance z(— as). The unit vector tangential to the 
axial helix is \a = dX/ds. The principal normal vector 
is fiA = (— cos(Kz), — sin(Kz), 0) and the binormal is 
i>A = tA x nA = a{sm(K z) , — cos(Kz), RK). Since tji is 
a unit vector sets a as: 



a 2 (R 2 K 2 + 1) = 1 



(1) 



The coordinates of the one-dimensional curve describ- 
ing the superhelix are: 

-R(s) = (R x ,R y ,R z ) 

— X(s) + r cos(ks)ftA + r sm(ks)bA ■ (2) 

Recall that the actual superhelices have nonzero thick- 
ness (the diameter of the copper wire), and hence are 
true three-dimensional structures. 



EXPERIMENT 

The classical experiments of Purcell, elaborated on in 
Q , examined the relationship between angular and trans- 
lational velocities of helical objects at very low Reynolds 
numbers. Here we extend these experiments to the su- 
perhelical objects described above. The experimental 
setup was originally designed for sedimentation experi- 
ments [HI (see Fig. QJa)). A tall transparent container 
is filled with silicone oil with large viscosity (y = 10 4 
cS, p — 0.98 g/cm 3 ). The oil behaves as a Newtonian 
fluid in the regime of interest here. Rather than allowing 
the superhelical object to descend by gravity, our exper- 
iment is designed to measure its rotational speed as it is 
towed up through the viscous column of fluid at a speci- 
fied translational speed. To drag the superhelix, a small 
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FIG. 2: Seven superhelices with increasing axial helix radius 
(from left to right). The upper panel is the side view and the 
lower panel is the axial view. 
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FIG. 3: A linear relationship between the translational veloc- 
ity and rotational frequency of a superhelix (r = 0.89 mm, 
R = 4.62 mm, p = 5.5 mm, and P = 19.4 mm). Triangles are 
from experimental observations. Dashed line is a least-square 
fit of experimental data. 



hook (~ 2 mm) is used to attach the superhelix to a 
thread from a motor. Note that the dimensions of this 
hook are quite small compared to the superhelix length 
(~ 4 cm) . By experimentally testing with an axisymmet- 
ric body (sphere), we found that this towing system (the 
thread plus the motor), does not produce any torque on 
the body. 

The superhelix is initially positioned near the bottom 
of the container, then is dragged upwards by the motor 
(Clifton Precision-North) at constant speed. In the in- 
termediate region in the container, steady state motion 
(constant translational velocity, rotational velocity, and 
drag force) is assumed. The superhelix positions, ori- 
entations, and velocities are measured from a 30 frames 
per second video stream of the camcorder. The trans- 
lational velocity in our experiments varies by changing 
power input to the motor. We have chosen a velocity 
range of 3-10 cm/s. Below 3 cm/s, the step motor pro- 
duces non-uniform pulsed axle rotations, which lead to 
irregular translational velocity. The Reynolds number 
based upon the towing velocity and radius of the super- 
helical structure (1 cm) is at most: 

Re= — ~0.1. (3) 
v 

We assume therefore that the steady Stokes equations 
govern the fluid mechanics of the translating superhelix. 
Within this translational velocity range, a linear rela- 
tionship between rotational velocity £1 and translational 
velocity U is observed (see Fig. [3]). 



A translating helix in a viscous solution rotates in the 
direction which it screws in. Following this rule, the small 
(straight) helix in our experiments would rotate clock- 
wise and the axial (straight) helix would rotate counter- 
clockwise when viewed from above. In Purcell's work [||, 
the jointed structure built by connecting two helices of 
opposite handedness, otherwise identical, showed no ro- 
tation during its sedimentation. The superhelix of inter- 
est here is the superposition of two helices with opposite 
handedness. The inherent rotational directions of these 
superimposed helices are in competition. For very small 
values of the non-dimensional parameter RK of the axial 
helix, the superhelical structure reverts to the straight 
small helix, and would rotate clockwise. One expects 
that for larger values of the parameter RK, the axial 
helix would be dominant, and the superhelical structure 
would rotate counter-clockwise. For some critical value 
of RK, we would expect a transition in direction, and 
hence, a structure that would show no rotation as it is 
towed through the fluid. We performed experiments that 
systematically varied RK, and observed this expected 
change in rotational direction. Fig. 2] shows the ratio of 
angular velocity to translational velocity as a function of 
RK, for the two different sets of superhelices (Sets I & 
II). Positive rotational rate is clockwise, and negative is 
counter-clockwise. In each set of experiments, the mea- 
sured ratio is depicted by triangles. In the next sections, 
we describe mathematical formulations that model these 
observations. 
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FIG. 4: Ratio of angular velocity to translational towing ve- 
locity. Triangles are values measured experimentally. Circles 
connected by lines are values predicted using the method of 
Regularized Stokeslets. Squares are values predicted using 
resistive force theory, (a) Superhelices of Set I. Positive rota- 
tional rate is clockwise and negative rate is counter-clockwise. 
Rotational direction changes around 0.7 for RK. (b) Super- 
helices for Set II. Same transition also occurs around 0.7 for 
RK. 



NUMERICAL RESULTS 
Regularized Stokeslets 

We assume that the superhelix is a rigid body moving 
in a Stokes fluid. The governing equations of motion are: 



- Vp + nV 2 u = 0, V • u = . 



(4) 



The total hydrodynamic force and torque exerted by the 
superhelix (with surface dD) on the surrounding fluid is: 



F = J f(x)dx, 



L = J (x — x ) x f (x)dx. 

xedD 



(5) 



(6) 



where f is the surface traction. 

A solution to the Stokes equations in 3D with a point 
force centered at Xq is the classical Stokeslet Due 
to the linearity of the Stokes equations, superposition 
of these fundamental solutions allows the construction 
of the velocity field induced by a distribution of point 
forces. The method of Regularized Stokeslets eases the 
evaluation of integrals with singular kernels by replacing 
the delta distribution of forces by a smooth, localized 
distribution [ll|, The force / = f(xo)5(x — Xq) is 
replaced by f — f(x )cj) e (x — x Q ) where e is a cutoff, 
or blob, function with integral one. This blob function is 
an approximation to the 3D Dirac delta function, with e 
a small parameter. Following [ll[, we choose 



,{X - Xq) = 



15e 4 



87r([[sc - SCq I 



,2\7/2 



(7) 



For N regularized point forces distributed on the sur- 
face of a body in rigid rotation and translation, the fluid 
velocity at any point x is evaluated as 



8nfiUi(x) =^2^2S- :j {x,x n )f j (x n ) 



For the given cutoff function, the kernel S is 



(8) 



2e 2 



(r 2 



e 2\3/2 



( r 2 + £ 2)3/2 



(9) 

where r = ||ac — x n \\. 

Note that evaluating equation ([5]) at each of the N 
points of the superhelix surface gives us a linear relation 
between the velocities and the forces exerted at these 
points. The matrix Sfj, for a given cutoff parameter e 
depends only upon the geometry of the superhelix. 

For a rigid body moving in a Stokes flow, there is a 
linear relationship between the total hydrodynamic force 
and torque and the translational and rotational velocity 
of the body 0. Following 0, [ll|, we focus on the z- 
components of total hydrodynamic force F and torque 
L, along with the z-component of translational velocity 
U, and rotational velocity about the z-axis f2. These are 
related by resistance (or propulsion) coefficients: 



= f 1 



A B 
B D 



(10) 



Here A, B, and D depend only upon the geometry of the 
object. 

In order to compute these coefficients, we discretize the 
cross sections of the copper wire using six azimuthal grid 
points. We choose a cutoff parameter e on the order of 
the distance between discrete points (see ll| for details) . 
At each point on this discretized superhelix, we impose 
a unit translational velocity and zero rotational velocity 
in the z-direction. We use the Regularized Stokeslet lin- 
ear relation (jHJl to solve for the forces on the superhelix 
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that produced this velocity. We then evaluate the inte- 
grals for total force F and total torque L in equations (0 
[5]) above. Using the linear equations (TlT)|) , we compute 
the resistance coefficients A, B. Similarly, we can com- 
pute D by imposing a unit rotational velocity and zero 
translational velocity 

These resistance coefficients allow us to predict the ra- 
tio of rotational velocity to translational velocity in our 
torque-free experiments described above, as 



O 

U 



B 
D 



(11) 



Figured] shows the Regularized Stokelet predictions (cir- 
cles) of these ratios for both sets of superhelices. The 
agreement with experimental data is excellent. Indeed, 
the transition from clockwise to counter-clockwise rota- 
tion is captured very precisely. 



Resistive Force Theory 



Resistive force theory [15j, [Iff is widely used to give 
an approximate description on a slender body moving in 
a viscous fluid. However, the non-local interactions of 
stress along the body are not taken into account. To see 
how important this non-local interaction is, we analyt- 
ically estimate the ratio of angular to translational ve- 
locities in this section, and compare the predictions with 
experiment and the numerical calculations using regular- 
ized Stokeslets. 

The Stokes drag force is proportional to its velocity as 
f = C t (u ■ t)t + C n (u ■ n)n where t and h are tangential 
and normal directions, and C% and C n are drag coeffi- 
cients. The normal direction is arbitrary, but, uniquely 
determined if the motion is given. 

First, we consider pure body-rotation about the z-axis. 
The body velocity is u = R±^ x tt = Q (R y , —R. x . 0) = 
ftR±u where R ± = (R x ,R y ,0), R± = \\R±\\, u = Rx 
z, and fi = Qz. The angle ip between the direction of 
motion and the tangential of the body is expressed by 
the superhelix coordinates as u ■ t — costp — (R y d s R x — 
R x d s R y ) j R± where t = R s . Similarly u ■ n = sin?/'. 

Using the vector relation ( ( A x B) ■ C = -(AxC)-B), 
the z-component torques associated with forces are 

L z (s) = (Rxf)-z 

= C t {u ■ i)(R xi)- z + C n (u ■ h)(R x h) ■ z 

= —C t (u-i)(Rxz)-i~C n (u-h)(Rxz)-h 

= -QR 2 ± (C t cos 2 tp + Cn sin 2 ip) , (12) 

and the total torque due to its rotational motion is 



Li 



Rotation) 



L z (s)ds 



-Q / R\ [C t cos 2 ijj + C n sin 2 ip] ds . 



Similarly, to total torque associated with the pure trans- 
lation u — (0, 0, U) — Uz is: 



ranslation 



> = U (C n - C t )(d s R z )R ± cosV-ds . (14) 



Decoupling the body motion into a pure rotation and a 
pure translation, the z-component of total torque on the 
body is expressed as 



L, 



■^(Translation) _|_ ^(Rotation) 



(15) 



In our experiments, we do not apply any external torque. 
Balancing two torques gives an expression for the ratio 
of angular velocity to translational velocity as 



-^(Translation) ^(Rotation) 

n _ J(C n -C t )(d s R z )R± cos ipds 
U ~ J R\ [C t cos 2 V + C n sin 2 V] ds 



(16) 



(13) 



Note that this ratio of velocities depends only upon the 
ratio of drag coefficients ^r-. Evaluating numerically 
these integrals, and using R and C„ = 2C t (the lead- 
ing order result of slender-body theory), Sl/U is plotted 
as squares in Fig. [4] 

The slope (-B/D) in the relation (fTT|) is measured in 
experiments by using a least-squares fit of the data (see 
Fig. [3]). Although the resistive force theory predictions 
of£l/U show the same general trend as the experimental 
measurements, the quantitative agreement is quite poor, 
and the transition from clockwise to counter-clockwise 
rotation is not captured. 



CONCLUSION 

In conclusion, we have studied the rotational dynam- 
ics of superhelices built out of superimposed helices of 
opposite handedness in a viscous fluid. As the radius 
of the axial helix increases, we have observed the tran- 
sition of rotational direction to the natural direction of 
rotation of the axial helix. The rotational velocities cor- 
responding to translational velocities are also predicted 
semi- analytically using resistive force theory, as well as 
numerically, using the method of Regularized Stokeslets 
flU fl2| . We see that although there is qualitative agree- 
ment between experiments and resistive force theory, for 
larger radii of the axial helix, the predicted ratios of ro- 
tation to translational velocities differ from experimen- 
tal results by more than one hundred percent, and show 
the incorrect direction of rotation. In contrast, the eas- 
ily implemented computational framework of Regularized 
Stokeslets demonstrates excellent quantitative agreement 
with experiments. 

While this study is motivated by the fascinating ge- 
ometry of spirochetes, we recognize that the rotation of 
the anterior superhelix of a spirochete is not due to an 
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imposed rotation about a vertical axis, but due to the 
counter-rotation of the periplasmic flagellum (that deter- 
mines the axial helix), against the cell body (the small 
helix). The fluid mechanic implications of this counter- 
rotation that govern spirochete motility is examined in 
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